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i—i ' Abstract 

> : 

, Rare radiative leptonic decay B — > Iv^ is studied in the instantaneous Bethe- 

Salpeter approach. The results are compared to other relativistic quark model 

■ predictions. 

o : 

The B-meson decay constant fg is an important phenomenological pa- 
rameter, not easy to measure directly. Purely leptonic decays B — > Wi, 
from which it could be extracted in principle, suffer either from helicity 
suppression m 2 /M| for light leptons, or from reconstruction difficulties for 

^ ; r-channel owing to the presence of two neutrinos in the final state. Decay 

• Eh ! rates expected in the Standard Model are 

X- 

"' nB^Wi) = ^\V ub ?fl,^Ml v L 3. 10 _ () . il = )l . 

where numbers quoted correspond to \V u b\ = 3 • 10~ 3 , fs = 200MeV and 
tb ~ 1.65ps. The present experimental limits Jlj on these decay rates are 
at least an order of magnitude larger. 

Some times ago Burdman, Goldman and Wyler (BGW) suggested an 
alternative, although model-dependent way for /# measurement [|2j|. The 
crucial observation was that the helicity suppression can be overcome and 
turned into an a (the e.m. fine-structure constant) suppression by an addi- 
tional photon emission in radiative weak decays B — > Wij. The dominant 
contribution in these radiative decays comes from the 5*-pole intermedi- 
ate state: a spin-0 B meson emits a hard photon and transforms into an 
off-shell spin-1 B* meson which by itself undergoes weak decay without 
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the helicity suppression. Thus measurement of these decay rates gives a 
tool to access B* meson decay constant f%. Heavy quark symmetry can 
be used then to relate and 

Afterwards B — >• Wij decay ( up to mf/Ml accuracy the decay rate 
is independent of the lepton flavor) was considered in a number of pub- 
lications. The BGW analysis was further refined by Colangelo, De Fazio 
and Nardulli j3j]. It was found, in particular, that the axial- vector (B') 
intermediate state contributes about 10% in the total decay rate. Atwood, 
Eilam and Soni used a simple nonrelativistic quark model to estimate the 
decay width and obtained || Br(B — ► Wij) ~ 3.5 ■ 10~ 6 , about 12 times 
larger than the purely leptonic branching ratio Br(B — > \xv^). Subsequent 
relativistic generalizations involve the use of: spinless Salpeter equation 
H, light front dynamics |3j|, light cone QCD sum rules [[7f, perturbative 
QCD combined with the heavy quark effective theory [|8| . All of them con- 
firm the main conclusion about the radiative decay mode enhancement, 
although with smaller branching ratio ranging from 0.9 • 10~ 6 to 2 ■ 10~ 6 
[j7j]. The recent experimental upper limits |9j are still far above of these 
predictions. 

It is worthwhile to note that the photon spectra in various relativistic 
quark models are significantly different as illustrated by Fig .[I]. An inverted- 
parabolic shape from ||, with a mean value ~ 1.3GeV, is asymmetrically 
modified in [|5], |5j towards higher photon energies, but in opposite direction 
in the light cone QCD sum rules approach [[7j , with the mean value shifted 
down to ~ 0.8GeV. This softening of the photon spectrum can effect the 
signal detection efficiency 0, [Tl| and the signal-background separation. We 
decided to check whether character of change in the photon energy distri- 
bution, predicted in 0, still persists in an another, instantaneous Bethe- 
Salpeter approach based, relativistic quark model QTT|, |T2j, |I3fl , which was 
successfully applied earlier to describe light |T|j and heavy |T3|, P~5| 



meson 



spectra, as well as electromagnetic form-factors of the ground state pseu- 
doscalar and vector mesons |TOJ , various two photon widths [[T7f and heavy 
meson weak decays [|TBj. This approach respects the heavy quark spin 
symmetry in the limit m& — ► oo [|T3j] - It also incorporates the spontaneous 



breaking of chiral symmetry in the light flavor sector, which happens due 
to the generation of a dynamical quark mass in the light quark self-energy 
generated by the interaction potential. 
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Figure 1: Photon spectra in various relativistic quark models for the B — > Iv^ decay. 
The B — > Wij matrix element has the following general structure 

A(B - w a ) = % ^^^ v "Wt^i - 7*>(p,) , 

where iif^, the hadronic tensor, is uniquely determined, due to gauge in- 
variance and Lorentz covariance, by two invariant form-factors Fy and Fa'. 

H, v = F A [g^P ■ k - k^P v \ + iF v ^ vaT P a k T (1) 

P being the 5-meson 4-momentum and k - the photon 4-momentum (in 
the following we will assume the 5-meson rest frame in all expressions, so 
F = (M B ,0)). 

The standard procedure leads to the following differential width (mi is 
neglected) 

d 2 T 
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B~ 



a{G F M B ) 
16tt 2 



\V ub \ 2 [p+( Xll x 2 )\F v + F A \ 2 + p_(x 1 ,x 2 )\F v - F A \ 2 ] (2) 



where p + (x\,X2) = (1 — 2x 1 )(l — 2x 2 ) 2 , p_(xi,X2) = (1 — 2a; 1 )(l — 2x 3 ) 2 and 
x\ = E^/Mb-, x 2 = E^/Mb, X3 = Eij Mb = 1 — x\ — x 2 are the photon, 
neutrino and charged lepton energy fractions. Fy and Fa form-factors 
depend only on the photon energy, so one integration can be readily done 
while calculating the decay width, and we get 



T(B -> Ivrj) = 



M B aiG ^ B)2 \V ub \ 2 j x\{l - 2 Xl ) [ \F v ( Xl )\ 2 + \F A ( Xl )\ 2 } dx r 



(3) 



The photon emission from the initial light quark gives the most impor- 
tant contribution to the B — > Ivfy decay amplitude 




k,e 

From this diagram the corresponding contribution to the hadronic tensor 
is easily obtained 



(4) 



-MbN c Q x 



dq 



2tt 



Sp {t m (1 " 75)G( 6 )(P + q)m P)G { u)(q)luG {u) (q + % 



where Q u = |, G^(p) = ^j- stands for the constituent g-quark prop- 
agator and T(q;P) - for the Bethe-Salpeter vertex function [|I2j| . The 
Mandelstam formalism [JTSj] or the bilocal effective meson theory [£2T| gives 
a general guidelines how to calculate transition or decay amplitudes, in- 
volving bound states, in terms of this vertex function. 

In the instantaneous approximation the Bethe-Salpeter vertex function 
V(q;P) depends only on the relative three-momentum q [|T2| . So the de- 
pendence on q in (4) is completely due to quark propagators and the 
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go-integration may be performed analytically using the residue theorem 
and 



G (q){p) 



Po ~ U( q ) (p) + it po + uj {q) {p) - 16 



7°, 



where u^(p) = yVn 2 + _p 2 and A±^ are the standard projection operators 
on positive and negative energies 



^±\P) = ^ (1± 
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After go-integration in (4) we obtain several terms, from which the lead- 
ing contribution comes from the following one 



H 



jlV 



dq M B N C Q U Sp | 7m (1 - 7s)r + -(g1 P)j„A^(q + % 



'27r) 3 [M b - u {h) (q) - cj (u) (q)] [k - u (u) (q) - u {u) (q + k)] 



, (5) 



because in the heavy 6-quark limit the characteristic momentum scale for 
r(<f; P) is much less than and so we may use 

M B ^m b ^ uj(b)(q) , U( u )(q + k) « /c > ^(u)® • 

In (5) T + _(g;P) = AV &) (g)7°r(g; P) 7 °A (u) (-0 and [t can be immedi- 
ately replaced by the Salpeter wave function <&(q) according to the Salpeter 
equation |T2|, 



r + _(g; P) = -\M B - co {b) (q) - ^©lA? 5 ®^®^^- 
In the heavy quark limit $(g) has the following simple form 

$(g) « -^L(l + 70)75^(5) , g = |g] , 



qV4n 



(6) 



where /(g) obeys the radial Salpeter equation and Srjj(q) is the inverse of 
the Foldy-Wouthuysen matrix for w-quark: 



1*1 



sin z/ (u) (g) 



the Foldy-Wouthuysen angle being determined as follows 



cos2z/ (u) (g) 



m. 



sin2i/( u )(g) 



^(u)(S 



To simplify (5), note also that 



A? } (£+ 2)7o = 5 ( "5(£ + q)\{l + 7o)5 (tt) (fc + g)7o = + q) + 7o)- 



But 



m, 



(A; + g) • 7 



^(«) + 0) u (u) (k + q) o; (u ) (k + q) 



and 



<%?(fc + 9)+7o 



1 



■[70 o;( u )(fc + g) - A; • 7] 



After these approximations (5) takes manifestly gauge invariant form: 



M B N C Q U f dq /(g) M7m(75 - 1)(1 - 70)^(4(07^} 



40F •/ (2tt) 3 q 



U( u ){q)u {u) {k + q) 



(7) 



Now it is straightforward to get Fy and Fa, the invariant form-factors 
from (7) as 

F v = F A = Q u (f(x 1 )+g{x 1 )) (8) 
where (we have introduced r = and x = dimensionless variables) 

oi V 2 

r + + # 



iV c IT 
= /o — / ^MI^O^ £i) - -^1)] 



(2tt) 5 /2 Xl 



p(xi) = 



00 



(2tt) 5 /2 3x? 



/(x)[^(x, — Xi) — ^(x, Xi)} 



\/r 2 + x 2 



Vr + ? — r 

vv* 2 + x 2 



<ix 



1/2 



— , (9) 



and 

</?(x,Xi) = 



\ 



r 2 + (x + X1) 2 
r 2 + x 2 



, ^(x, Xi) = (r 2 + x 2 + x 2 — xxi)(/?(xi) . (10) 



Note that in the nonrelativistic (heavy w-quark, but m u <C ra&) limit 
characteristic momentum scale for /(x) is much less than r. So g(xi) — > 

—* ^2^)5/2 ^7 j xl{x)dx . 
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In this limit 5-meson decay constant fs takes the form \\T3 



30 



B 



So 



» 5 A 



Mb \ xl(x)dx 



1 / 



B 



2x\ m u 

Substituting this into (8) and (3), we reproduce Atwood, Eilam and Soni's 
result 



a(G F M%f 2 2 f% 

•y u \vvb\ 



288tt 2 ' 

The equality of the vector and axial current form factors, given by 
is surprising because in the pole approximation they receive contributions 
from intermediate states with opposite parities. Nevertheless this addi- 
tional interesting spin symmetry was shown to hold in the high recoil 
E 1 ^> Aqcd region at least at one-loop order by explicit leading twist 
perturbative QCD calculation ||. 

The radial wave function l(x) is determined from the integral equation 
which in the heavy 6-quark limit takes the form |T3 



[M -m b - u {u) (p)} l{jp) = -J dq 



\ 



1 + 



w (u)(p). 



1 + 



\ 



1 



ra. 



1 



m 7 



•vi(p,q) 



Kq), 



>v (p, q) + 



where M is the bound state mass and vi(p,q) angular matrix element of 
the potential kernel is determined through 

PI 
(2ny 



dVL p / d£l q Yl, M ,{p)V{p- q)Y LM {q) = v L (p,q)5 LL >5 M M> 



For estimation purposes we have used the same parameter set for the B- 
meson description as given in jT3j . That is constituent quark masses m& 
1 . 79 GeV, m u = 0.33 GeV and the linear plus Coulomb potential 

V(r) = --— + o 2 r, 
6 r 

with a s = 0.39 and a = 0.41 GeV. 

For these parameters the radial Salpeter equation flITf) was solved by 
Multhopp method [p~3| , |2Tf . The resulting radial wave function is shown in 
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Figure 2: The radial wave function l(x). 

Fig.0. Note that ( [TT|) determines l(x) only up to normalization constant. 
The normalization condition was considered for a general case in |T3j and 
for our approximations looks like 




l 2 {x)dx = 1. 



Having at hand the radial wave function, we can calculate the decay 
branching ratio 

Br{B li>a) « 0.9 x 1(T 6 . 

The resulting photon spectrum is indicated in Fig.rTJ- As one can see, our 
results are very close to predictions of the spinless Salpeter equation model 
j5j , except in the low-energy part of the photon spectrum where the results 
of j5j] contain unphysical divergence and are not reliable. 
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To summarize, the instantaneous Bethe-Salpeter model, considered in 
this article, gives a photon spectrum similar to other relativistic quark 
models |5j and j6j], but different from the light cone QCD sum rules ap- 
proach j7j] . The predicted branching ratio is within the reach of the BaBar 
experiment So we may expect that some experimental information 

will appear about this rare decay mode in the near future. To extract inter- 
esting quantities like fs from this information, a detailed understanding of 
the model uncertainties in simulation of this decay is necessary. We hope 
that our investigation will be useful in such studies. 
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